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We study the O(N) quantum Heisenberg antiferromagnet using a parametrisation in terms of real 
fcrmions. The N — > 00 limit of the unfrustrated model is controlled by a saddle point representing 
singlets on dimers that cover the lattice and which is infinitely degenerate in many cases. An 
infinite degeneracy of the N = 00 ground state occurs also in some, but not all, frustrated quantum 
antiferromagnets. Our results are similar to results reported previously on an SU(N) generalisation 
of the Heisenberg antiferromagnet for N = 00. 
PACS numbers: 75.10.Jm, 05.30.-d 

INTRODUCTION 

The idea that the ground state of a frustrated spin 1/2 antiferromagnet should be describable in terms of products 
of spin singlets was suggested originally in the context of triangular spin | antiferromagnets Q| . Products of singlet 
states are exact ground states of certain frustrated spin 1 /2 antiferromagnetic spin chains with special couplings 0] 
and they have been used in variational approaches |3j and in heuristic arguments. 

In the present note we study the O(N) Heisenberg antiferromagnet in the N — > 00 limit by using a constraint free 
parametrisation of the O(N) Heisenberg antiferromagnet in terms of real fermions. We show that the N — > 00 limit 
is dominated by a mean field configuration with small ~ fluctuations. We are able to locate the saddle point of 
lowest energy in the case of unfrustrated antiferromagnets and for a subclass of frustrated ones. However, in many 
cases we do not find a single saddle point, but instead an infinity of saddle points that are all degenerate in energy. 
These saddle points are essentially singlets on nearest neighbor dimers that cover the lattice and correspond to the 
RVB singlets of Anderson. 

Our results are similar to those obtained by D.S. Rokhsar on the N — > 00 limit of the SU(N) extension of the 
Heisenberg antiferromagnet M. 

MAJORANA REPRESENTATION OF THE O(N) HEISENBERG MODEL. 

It is known that operators that operate on spinors can be generalised from the 0(3) group to O(N) by using O(N) 
Dirac matrices. Somewhat less familiar is the fact that these Dirac matrices can also be considered to be real fermions 
Hi, thus providing us with the following representation of O(N) spin matrices: 

1 

Sk.i = - ^ (VkVi ~ VlVk) ,k,l=l,2, ..N (1) 

[vk,vi]+ = s ki> ? // + = m 

The representations of Sk,i can be characterised by the value of their quadratic invariant 

\ Yl s v s v=-\ ViVjViVj = I (N 2 ~ N ) (2) 



s 2 

i,j=l..N i^j=l..N 



For N = 3 and 0(3) this representation reduces to 

Si = -irtmu i,k,l = 1, 2, 3 cyclic (3) 
9 1/1 
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a representation that has been discussed extensively in the first of ref 0. To define the O(N) Heisenberg model, we 
associate a spin represented by fermions with each point of a lattice: 



S x ,kl = - ^ (VxkVyl - Vyl^xk) 
[VxkiVyl] = $klSxy 



(4) 



and use the invariant scalar product 



a? 



(SxSy)o(N) — 7} J]] s x,kl ■ Sy^kl 

to define the interaction: 

H H eisenber g — ^ ^ Jxy^x^y)o{N) ^ ^ ^ 



(5) 



<x,y> 



<x,y> 



N 1 
J + 2 



I N \ 



\k=l 



(0) 



The fermions rj x ± are real as is appropriate in a situation without electric charges, and there are no constraints on the 
physical Hilbert space. To get familiar with the representation of spins in terms of real fermions, we may calculate 
the energy of a pair of nearest neighbor points x, y or a "dimer " that interact via their spins. We do this most easily 
by introducing ordinary complex fermions: 



^-dimcr,fc — (Vx,k + i^y,k) , ^dimcr.fc 



1 



V2 



(Vx,k - iih.k) , k = 1, ..N 



(7) 



(SxSy) 



N 1 



O(N) ~ g 



A' 



(^dimcr , k ^dimer , k ^dimcr , k ^ dimer , fc ) 



,fe=l 



and find 



N 2 N 



^dimer — < ( s xS y ) ^ N ^ >— + 

The ground state energy of a dimer is seen to vary smoothly as a function of N. 



(8) 



A TOY MODEL 

To gain further insight into the nature of the N — > oo extrapolation, we calculate, for arbitrary N, the ground state 
energy of four spins located at the corners of a square with interactions along the edges and across the diagonals: 

H square = SlS 2 + S 2 S 3 + S3S4 + S4S1 + £ (S1S3 + S 2 S 4 ) 

= \ [(«1 + S3) + («2 + (9) 



+ £ ~^- ((Si + S3) 2 + (s 2 + S 4 ) 2 ) 



2es 2 



Since one knows how to combine representations of 0(3), one easily obtains the spectrum of H square for N — 3. One 
finds a singlet for the ground state with an energy of 

E S quare,3 = j £ 3e \ £ ^ | (10) 

We assume in our calculation for O(N), with arbitrary N > 3, that the ground state continues to be a singlet, so 
[(si + S3) + (S2 + S4)] 2 = in eq(^|). It remains to find the spectra of (si + S3) 2 and (s 2 + S4) 2 . We now use standard 
Dirac Gamma Matrices instead of real fermions in our argument and reduce the tensor product of two O(N) spinor 
representations, say ^ a (fP , of 2^ N ^ 2 ^ dimensions each, to calculate (si + S3) 2 and (s 2 + S4) 2 . The reduction of £ Q <^ is 
achieved via antisymmetrised products of Dirac O(N) gamma matrices 
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eV - ft, er M1Pa *7, ...,ar^. m r, (11) 

where r„j ,„ n denotes a product of n gamma matrices that are antisymmetrised in their indices. The O(N) spin 
of an antisymmetric tensor r^..^ is most easily found by treating it as Grassmann numbers acted upon by O(N) 
generators = Tj — Tj -^jk- . In this way one can find without too much difficulty that 

s 2 ^..^ = i ^ s ij s ij r iil .. lin =n{N-n)T lil .. IMn (12) 

i,j=l..N 

with a minimum eigenvalue of s 2 at n = [^] = T — 4^iV,odd an d with maxima at n = 0, AT. This results in the following 
ground state energy of H square for arbitrary N: 

^are,iV - | _ 1^2 + 1^ ^ j j (13) 

The ground state enery E square ^ in eq(|l3|) reduces smoothly to its value at A" = 3, with a relative precision of order 
-^2 ~ 10% for N = 3, if one includes only the two leading terms, and we conclude that the 0(N) extrapolation is 
satisfactory in this toy example. By comparing with eq(|^) we notice that the leading 0(N 2 ) term of E square ^ can 
be interpreted in terms of the formation of two dimers with coupling s or 1, whichever is larger. We shall see below 
that this is part of a more general pattern that emerges at A" = oo. 



NATURE OF THE SADDLE POINT AT N=oo 
To use standard N — > oo techniques || we rewrite the quartic interaction in terms of an auxiliary field 
H = Y J ^0* s 2/-y)= Y '-^-{VxVyf 

<x,y> <x,y> 

Drje' So dt (i^9 t ^+H) = cmst{l3) J DBDrje -S 



const(f3) 



J DBe So dt{ * E.., ^"f u»g(det(a r +ifl)) 



f f 1 -D 1 % \ 

s = dt ( 1 Y + 2 Y ^vPMw + 2 X! />> < •/'/.;. '/,//• j (14) 

With scaling as the exponent is of order A/" and the integration over is dominated by "classical " config- 
urations of B xy plus fluctuations of B xy of order 0(1/ V~N). To find the optimal solution(s) of eq(pT[) at T = it is 
easiest to consider the Hamiltonian that describes the saddlepoint of eq(|14[) and which is given by 

H °° = I Y "T 1 + \ Y iB *v r lxn'nvi* (1 5 ) 

x,y y x 7 y,iJ, 

The optimal auxiliary fields B xy are those that give the lowest energy and which minimise < >. The matrix 
iB x ^y is hermitean and antisymmetric and its eigenvalues occur in complex conjugate pairs (£ n , A n ), (£*, — A„), for 
n = 1, 2..\H=(points). This property of iB XtV enables us to rewrite the kinetic energy of the real fermions r/ xfl in terms 
of complex fermions with standard oscillator commutators: 

i#(pomts) at 

2 iBxyTlx^rjy^ = — Y^ ^™ (^n,(U^n,A« — (1^) 

x 7 y,fi n—1 fi— 1 

With i4„ jM = ^ tnxVxfi, 
x=l 

A ni „] + = <W<W; [A mi p ) A n , v ] + = 0]m,n = l..~#(points) 
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We deduce from eq(16) that the zero point energy of the fermions is given by 



< \ iB xyVx^Vyn >= l A l 

x,y,\i all A 

where one must sum over all ^(points) eigenvalues. Combining eqs( jl5|Jl7| ) we may rewrite the ground state energy 
at N — 00 as 

X r / Tit j: 

x,y Xtf all A 

In the case of frustrated magnetism we have, in general, distinct matrices and iB that do not commute. However, 
for a single nonzero coupling J xy — J the energy associated with the auxiliary field can be simplified and related to 
the spectrum of the matrix iB: 

53 "T^ = 7 ^2 B *v = — j Y B xy B yx, B yx = -B xy (19) 



x,y 

.-TrB 2 = -Tr (iB) (iB) = - V A 2 

J J aUA 



Combining cqs(|15|,|l7|,|l9|) we obtain a lower bound on the ground state energy: 



V /A 2 \X\N\ 1 ^ 

I 4 7 4 / 4 7 Zv 



4J 4 / 4J 

all A v 7 all A 



2 (J7V) 2 



21 7 \r2 



JiV 2 

> -^-#(points) (20) 



We have used that the matrix B x y has as many eigenvalues as there are points on the lattice. We have thus obtained 
a lower bound for the energy of an O(oo) antiferromagnet with a single coupling J. 

For unfrustrated Heisenberg Hamiltonians on a lattice of points that can be considered as a union of non overlapping 
pairs of points or "dimers " it is easy to saturate this bound. For such Hamiltonians, we may choose a configuration 
of B xy that is nonvanishing only on an arbitrary collection of nonoverlapping dimers that cover the whole lattice. By 
hypothesis, the couplings are the same on all these dimers. The matrix B XiV then decomposes into blocks, one for 
each dimer, and according to eq(l2C|) the energy can be rewritten as 

2 ' 4" 



jt 1 V" \(\\\ JN \z (JiV) 2 ] |A|=^ 1 ^ ^ (JN) 



4J 

dimers 1—1.2 

JN 2 „ JN 2 



4J 

dimers 2—1,2 



• 2 • #(dimcrs) = — • #(points) (21) 

16 16 

Here we have adjusted the block of B that corresponds to each dimer in such a way that its eigenvalues are A = ±4^-- 
We may also return to equation (||) to see more directly that the minimal energy of a collection of nonoverlapping 
dimers (in their singlet state) saturates the inequality (|20|): 



JN 2 

J(sxSy)o(N) >= g- ■ #(dimers) + O(N) (22) 

' x,y> 



< 

<x,y> 

which visibly saturates eq(pl]). Returning to the toy model of eq(j^) we now understand why its ground state energy 
corresponds, to leading order in N, to that of two dimers, but our general arguments apply only to e — 1 and e — 
where there is a single coupling. The toy model suggests the stronger statement that the ground state on lattices 
coverable by dimers is one of singlets on the "strongest " dimers. 

To follow the hint of the toy model, we reconsider eq([l8]) and derive a lower bound for the ground state energy of 
a frustrated antiferromagnets at N = oo: 



4 x,y Jx V 4 all A 4Jmax x,y 4 all A 

Here we have used the positivity of the couplings J xy . We may now copy word for word the arguments that lead to 
the lower bound of eq(pO|) in the unfrustated case, but with J replaced by J max - We find: 
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(24) 



In some cases this bound can be saturated and lowest energy saddle points can be found. Consider, for example, 
the square lattice spin 1/2 antifcrromagnet, with nearest neighbor couplings J\ and next nearest neighbor couplings 
(across the diagonals) J 2- In this case we chose the stronger coupling J max = max(Ji, J2), and any dimer covering 
of the lattice by the "stronger bonds " that correspond to J max . As before, we saturate the lower bound with these 
configurations. So the infinite degeneracy persists in this particular frustrated antifcrromagnet. The N = 00 saddle 
point of an antifcrromagnet on a Kagome lattice Q is also infinitely degenerate, because a Kagome lattice can be 
covered by dimers in an infinite number of ways. However, there is only finite degeneracy in the N = 00 saddle point 
of the frustrated spin chains of (|] . 

Although we have found an infinite number of degenerate saddle points in certain O(oo) antiferromagnets we cannot 
be sure to have found all the saddle points that saturate the lower bound. In particular, the saddle point may exhibit 
continuous degeneracies. A one parameter continuous degeneracy is indeed present in the O(oo) saddle points of the 
frustrated square of eq(^|), and an analogous degeneracy was also noted in || for the case of SU (00). 



CONCLUSIONS AND OPEN PROBLEMS 



We have found that a large class of O(N) Heisenberg models have N = 00 saddle points consisting of products of 
singlets and which are infinitely degenerate in many cases. 

Because there is a gap at large N but no gap in the unfrustrated model at N = 3 and N = 4 there should be a 
critical value of N that separates the two regimes in the case of unfrustrated antiferromagnets, while for frustrated 
antiferromagnets there is no need for such a phase transition. 

Our results closely parallel those obtained by D.S. Rokhsar [0 on the N — > 00 limit of an SU(N) extension of the 
Heisenberg antifcrromagnet. In the latter case, the effects of ljN corrections are understood 0, while in the O(N) 
case, the effect of these corrections is still unknown. Also, it is still an open problem whether the correlations can be 
usefully organised in powers of N at N = 3. 
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